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At  the  present  time,  semiconductor  thermobatteries  are  used  as 
thermostats  in  radioelectronic  d  vices,  air  conditioning,  refrigerating 
of  food  products,  biological  preparations  and  the  like.  The  operation 
of  thermoelectric  devices  is  established  by  both  static  characteristics 
(refrigerating  capacity,  fixed  temperature  drop)  and  by  dynamic  ones 
(exit  time  from  a  given  temperature  mode,  precision  of  maintenance  of  a 
temperature).  Of  special  interest  among  possible  transitional  modes  of 
operation  of  thermobatteries  are  transitional  processes  arising  during 
passage  of  multistep  currents  across  thermoelement  junctions.  A  number 
of  variants  of  calculations  of  temperature  movement  when  thermobattery 
power  is  switched  on  have  been  published  in  the  literature  [1-8].  As  a 
rule,  simplified  thermobattery  models  are  discussed  in  the  works  cited. 
Either  they  generally  do  not  take  the  therraoload  on  thermobattery  junc¬ 
tions  into  account  [3,  4],  or  they  only  take  into  account  the  heat  capa¬ 
city  of  the  switchp.late  and  the  object  cooled  [1,  2,  5]. 

Calculation  of  transitional  thermal  processes  on  a  thermoelement 
cold  junction  under  mixed  thermoloading  is  presented  below.  Subsequently, 
it  is  considered  that  thermobattery  junctions  have  identical 

electrical  and  thermal  characteristics,  and  that  the  side  faces  of 
thermobatteries  are  adiabatically  isolated.  Then,  temperature  distribu¬ 
tion  T  of  a  thermobattery,  as  a  function  of  coordinate  x,  directed  along 
the  length  of  thermoelements  from  the  cold  junction  to  the  hot,  and 
time  t  are  written  in  a  thermoconductivity  equation  as 


Here,  X  is  the  thermoconductivity  of  thermoelement  junctions,  c  is 
the  volumetric  heat  capacity  of  thermobatteries ,  p  is  the  resistivity  of 
thermoelemen*  junctions,  j  is  the  current  supply  density  (here  and  further, 
the  values  of  physical  parameters  are  related  to  unit  areas  of  the 
thermobattery  cooling  surface) . 


(2a) 


!  K=„  “  •'0  + ^TP57  +  B<.  -  »o)  +  /<„ . 

«|(„=V  <3a) 


Using  the  operator  method,  it  is  easy  to  find  the  Laplace  transform 
for  the  temperature  of  the  cold  junction 


n\\  -  °o  .  'V/  (ch  VP  —  l)  —  p  (X.,  —  a:,,.)  sJi 

(P)U=0  P  +  /  \\SJ  Ch  YJ  4-  (v,B!’;-W)  sh  YJ\  ‘ 


(4) 


Transition  to  the  provisional  domain  is  accomplished  by  the  Rieman- 
Mellin  formula  and  residue  theoiy: 


..0  0,5 ,2-  Kr. 

0(Fo)  =  0(1  —  . i  TTTW 


^  -4,exp(—  5;fo). 


(5) 


where 


^  _  o  1  (  -t>0  —  K„)  S'''')  ?v  —  ( 1  —  cos  Vl  I _ 

5k15k(2,'i  “  1 ) s,n  I 5k—  1  -  '  -  Bij)  cos  0,. 


The  value  of  6K  is  the  positive  root  of  the  transcendental  equation 

tg  0  “  T  iO- - V  -  Ul  i  ’  (6) 


It  follows  from  formula  (5)  that  the  temperature  drop  A9=0q-6  between 
thermoelement  junctions  can  be  presented  in  the  form 


where 


A0=,A0CI,«-/(Fo). 


/(Fo)=-  2  A.cxp(-o‘Fo). 

k  r=  i 


(7) 
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Series  (7)  defines  the  rate  of  approach  of  the  temperature  of  the  cold 
junction  to  a  stable  value,  in  which 


AOc»u  =■ 


•'8u- (0.5 

1  +  v  +  Bl] 


(8) 


The  maximum  temperature  drop  is  achieved  at  a  current  density  v®,  the 
value  of  which  is  equal  to 


o  \B-  -!-  20o/J  H-  2/Cj,  (1  -}-  2£K)l  ^  -£ 
Vr~:  1 


(9) 


in  which  the  value  itself  of  the  maximum  drop  is  determined  by  the 
following  relationship: 

A'£au=  0„ (10) 


here,  B=(l+Bii) (1+2CK) . 

Calculations  according  to  formula  (5)  assume  knowledge  of  the  solutions 
of  Eq.  (6).  The  values  of  the  two  least  roots  6j  and  62  were  calculated 
on  a  digital  computer  for  a  wide  range  of  changes  in  the  parameters 
v+Bij  and  m-  The  results  of  the  calculations  are  presented  in  Table  1. 

The  case  for  v+Bij-0  and  m=0  are  limiting.  Since  the  value  of  is 

localized  in  the  corresponding  intervals  6K£((k-l)7r,  kn)  ,  series  (7) 

quickly  converges  at  sufficiently  large  values  of  Fo. 

It  is  evident  from  the  data  of  Table  1  that  the  increase  in  heat 
capacity  m  leads  to  a  decrease  in  the  values  of  6  ,  that  is,  to  an  increase 

in  retentivity  of  thermoelements;  increase  in  heat  transfer  Bij,  as  well 
as  current  density  \>,  show  an  opposite  effect  on  the  value  of  6^. 


Pre-exponential  factor  A^,  also  depends  on  parameters  m: 
In  addition,  the  parameters  0o  and  influence  their  values. 

9  2 

of  changes  in  coefficients  Aj  and  A2,  as  well  as  of  6J  and  S2 


v  and  Bij . 
The  course 
vs.  the 


ratio  of  current  density  v  to  its  optimum  value  v®  for  various  combinations 

of  parameters  Bi*  and  m>  are  presented  in  Fig.  1,  2,  3  and  4.  The  rela¬ 
tionship  of  Aj  to  current  v  is  characterized  by  the  presence  of  a  maximum, 
after  which  Aj  decreases  monotonically ,  passes  through  zero  and  becomes 
negative.  Coefficient  A2  increases  monotonically  with  increase  in  current 
v.  With  increase  in  heat  transfer  in  Bij,  the  value  of  the  coefficients 
of  the  series  decreases,  which  corresponds  to  a  decrease  in  the  steady 
drop  with  increase  of  Bij  (Fig.  1  and  2).  Values  of  and  52  increase 
monotonically  together  with  current  strength. 
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efficient  Aj  vs.  current  density  v 
i0=0.6) 

Bij=0,  r,!  =  l 
Bii=0,  ni=0 
Bii=5,  m=0.1 
3ii=5,  m=0 
Bii=5,  m  =  l 


In  those  cases  when  the  parameters  Bi;  and  m  assume  limiting  values, 

close  to  zero  or  fairiy  large,  it  is  easy  to  determine  analytically  the 

corresponding  values  of  6  and  A„. 

K,  K. 

Thus,  when  v  -f  Bh-*-0,  K»=0 

when  v  +  Bii»l,  K„= 0 


;  - k -  A  -  2v2|(-l)*-H 

K  «.L,  /1K—  -2;.2  (v  .1.  Bi,)  ’ 


ni»l.  Kn—Q 


=  (*  —  1 ) ~  (k=-2,  3 . ), 

4  -_:4-0,5^  _.2vall  +  (-l)*i  (k  —  0  3  ) 


Correspondingly, temperature  changes  at  sufficiently  high  values  of 
Fo  are  determined  by  the  following  relationships: 

when  n, -3- 0,  v  +  Bii-^0,  /C„=0 


A0~*  AO, 


when  V  +  Bii»l,  /C„=0 


A0  =  AOc»u  +  ~ (T:7'Bii)  exp  (---  Fo), 


when  r)i3>  1 


AO  AGctju [l  -exp(--i±^->-}i|  Fo)]. 


The  temperature  of  the  cold  junction  in  the  transition  mode  of 
operation,  depending  on  the  current  strength,  can  approach  equilibrium 
from  the  sides  of  large  or  small  values,  that  is,  monotenically  or 
passing  through  a  minimum.  The  condition  fixing  the  limit  for  these  two 
cases  is  the  equality  of  coefficient  Aj  with  zero: 


•-(».-  *•) 


sino,  —  v(l  —  roso,)  -0. 


Determination  of  the  current  vq  which  satisfies  condition  (14)  should 
be  accomplished  numerically  together  with  the  solution  of  Eq.  (6).  It  is 
interesting  to  compare  the  values  of  vq  obtained  here  with  the  optimum 


. . . . . . . . . . . . . . . tel®:*  . . . . . . . . . 


/ttfr-m-Vi'l  f  A' 


ones  in  the  steady-state  mode  of  operation.  Data  of  the  calculation  of 
vq  for  a  series  of  values  of  the  parameters  m  and  Bij  (at  0q=O.6  and 
K^=0)  are  presented  in  Table  2.  When  -*■  0,  the  value  of  Vq  can  be 

found  analytically, 


60arccos 


Jli— ' \ 

3  I  1  -r  «o  f 


V  1  -i-2U„ 


Since 


v°t= y  l+20o  —  1 ,  then 

=  -*-  fl  -j — 7====-)  arccos  (—  -r g-) . 
.0  2  ’v  ^  1^1  2»o  I  \  1  +  b0  j 


At  any  0q,  the  latter  ratio  exceeds  -j,  that  is,  it  is  greater  than 

unity.  An  increase  in  the  thermocapacitive  load  m  favors  an  increase  in 
current  v0,  but  heat  transfer  Bi*  decreases  the  value  of  vq  (Table  2). 

v0 

In  the  general  case,  the  ratio  — ij-  can  be  both  larger  and  smaller  than 

VT 

unity.  Thereby,  the  extreme  on  the  curve  of  the  function  0(Fo),  arising 
at  v=vq»  can  appear  at  current  values  lower  than  the  optimum  (at  fairly 
large  values  of  Bi])  and  in  excess  of  the  optimum  in  the  steady-state 
mode  of  operation.  Since  62:>6i,  the  mode  of  operation  in  transition  from 
a  monotonic  temperature  drop  with  time  to  a  temperature  change  with  the 
extreme,  is  the  most  "rapid."  Current  density  v0  corresponds  to  the  shortest 
time  for  the  transition  process,  and  the,v=VQ  mode  of  operation  is  advisable 
in  conditions  when  the  current  vq  provides  a  sufficient  lowering  of  the 
temperature  in  the  steady-state  condition. 


Table  2:  Relationship  between  current  vq  and  at  various  thermal 
loads  on  the  cold  junctions 


T.l  - 

... 

.  0 

0,1 

0.5 

1 

oo  0 

0 

0 

Bl, 

.  .  . 

.  0 

0 

0 

0 

0  0 

1 

5 

0 

VT 

... 

.  0.-1&3 

0,483 

0,483 

0,483 

0.4S3  0,433 

0,530 

0 ,373 

^0  • 

.  .  • 

.  0,79 

0.83 

0,97 

1,05 

1 ,20  0,79 

0,« '6 

0.37 

V  = 

V?  • 

.  1 ,64 

i.74 

2,00 

2,17 

2,4S  1,61 

1,25 

0,65 

r.i  ■ 

.  0 

0 

0 

0,1 

0,5 

1 

A, 

Bl, 

.  .  . 

.  10 

cc 

5 

5 

5 

5 

5 

.  .  . 

.  0,584 

0,600 

0,573 

0,573  0,57.3 

0,573 

0.573 

'■'a  • 

.  .  . 

.  0,24 

0 

0,37 

0,41 

0,55 

0,72 

1 .20 

v  = 

•'0 

-  0,41 

0 

0,65 

0,71 

0,97 

1,25 

2,09 

In  solving  the  termal  conductivity  Eq.  (la)  in  form  (5),  a  rapid 
convergence  is  observed  for  sufficiently  high  temporal  values  of  Fo. 

For  plotting  the  transition  process  at  small  temporal  values  of  Fo, 
the  solution  of  Eq.  (la)  should  be  used  in  the  form  of  an  exponential  series 
Fo.  From  formula  (4),  the  expansion  of  0(p)  by  exponents  p-1'2  can 
be  obtained: 


hi  „i  %  ‘‘8|| — 1  i  v«2  +  '-Oo  —  Kn 

Hp)--p - —  ~^  + - x - 

(•<Cp  —  Kh)  II  —  T,i  (•'  t  Bh)  I  ~i~  T-iv2  | 

r.3  ' 

'• \P 

Wo—K„)  II  -  2m  ('4  Bi,)l  +  M'2  H  —  T.i  i-  Bl))l 


which,  in  the  temporal  domain,  corresponds  to  the  series 


0 (Fo) l-o 


3  V  ~- 

(v«0 ~  ATh)  1 1  -r„  (v  -r  Bl,)l  + 


T0 


Fo2- 


Fo2  - 


■X 


2-?  1  i5i^r 

x  (  |1  -2r,,  (y-h  Bi,)|  -j-  r.rj2  1 1  —  T.i  (v-f  Bl,)l  Fot 


(16) 


Similarly,  the  expansion  of  6(Fo)  by  powers  of  Fo  can  be  found  for  the 
case  of  a  negligible  heat  capacity  qj : 


0 (Fo)  —  O0  -  -pL-(v0o  -  !<„)  V  Fo  ■«-  [(*0o-  /C„)(v  +  Bi,)  + 

4-  v2]  Fo  -  -^=-  [(V0O  -  K„)(y  ~r  Bi,)  +  v2]  (v  +  Bi,) FoT  + 
4-  4-  IC4  -  *  „) (v +  Bi,)  +  v2I  (v  +  Bi,)2  Fo2  -  . . . 


(17) 


Comparison  of  formulas  (16)  and  (17)  shows  that  the  presence  of  a  thermal 
capacitance  on  the  cold  junction  leads  to  a  slowing  down  of  the  rate  of  cool¬ 
ing  in  the  initial  period.  The  first  two  terms  of  expansion  (17)  coincide 
with  the  corresponding  terms  of  the  expansion  for  the  case  of  cooling  of 
the  plate  by  a  constant  output  source,  equal  to  vep-K^  [9]. 

To  evaluate  the  accuracy  of  calculations  of  temperature  by  the  analytic 
formulas,  using  the  first  terms  of  the  corresponding  series,  the  data  of  the 
numerical  solution  of  Eq.  (3a)  can  be  used.  A  graph  of  change  in  temperature 
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drop  between  battery  junctions  with  time  for  various  values  of  the  ratio 
an,1  the  parameters  m  and  Bij  (6o=0.6,  K^=0)  is  presented  in  Fig. 


v  = 


In  the  case  m=0,  Bij=0  (Fig.  5a),  large  current  values  corrsspond  to  the 
curves  with  large  angles  of  ascent  at  the  coordinate  origin.  Curves  for 
currents  less  than  vq  (for  the  case  under  consideration,  vq=1.64v^,)  have  a 

monotonic  character  and  do  not  intersect  when  v<v®,  but.  curves  corresponding 

to  currents  exceeding  the  value  of  v®.  begin  to  intersect  the  preceding 

ones.  When  the  current  exceeds  the  value  v0,  the  temperature  trend  has 
an  extreme.  Curve  5  (v=1.5v°  and  near  to  the  value  of  vq)  has  the  greatest 

rate  of  establishment  of  the  steady-state  fall  for  the  relationships 
presented. 

The  graph  in  Fig.  5b  (m=l,  Bii=0)  makes  it  possible  to  follow  the 
effect  of  thermocapacitive  loading  on  the  course  of  cooling.  An  extreme 
is  not  observed  in  the  curve  A0=f(Fo),  even  at  v  2  (for  this  case, 
vo=2.17v®.  The  fastest  rate  of  establishing  the  steady-state  condition 

among  the  curves  presented  corresponds  to  7-2. 

For  greater  values  of  the  heat  transfer  coefficient  (Fig.  5c,  Bii=5, 
*11=0),  the  value  of  7  at  which  an  extreme  temperature  trend  appears  is  less 
than  unity,  and  the  curve  A0 (Fo)  for  v=v°,  has  an  extreme.  The  temperature 

drop  at  the  extreme  points  of  the  curve  A0(Fo)  exceeds  the  steady-state  one. 

The  case  of  Bi*=5  and  m=l  (Fig.  5d)  is  characterized  by  mutual 
compensation  of  the  effects  of  heat  transfer  and  thermocapacitive  load, 
as  a  result  of  which  the  extreme  on  the  curves  arises  at  7=1.25 

The  data  from  calculation  of  function  A0(Fo)  according  to  formula  (5) 
are  noted  in  circles  in  Fig.  5a,  b,  and  c,  and  by  formulas  (16)  or  (17), 
by  triangles.  A  comparison  with  curves  plotted  according  to  the  results 
of  calculation  of  function  A0(Fo)  by  numerical  methods  on  a  digital  compu¬ 
ter  shows  that,  for  Fo£(0;  0.05),  good  precision  of  the  approximation  is 
insured  by  three  terms  of  series  (8)  or  (9),  and  at  Fo>0.04,  sufficient 
accuracy  is  achieved  with  the  calculation  of  two  terms  of  series  (5). 

The  formulas  obtained  above  for  calculation  of  the  cooling  processes 
can  be  used  in  discussion  of  the  transition  processes  in  the  case  of 
heating  (if  the  temperature  of  the  cold  junction  is  held  constant  here), 
substituting  -v  for  the  value  of  v  in  the  corresponding  formulas.  The 
parameters  correspondingly,  will  be  the  roots  of  the  equation 


t"  3: 


,4- 


Dii 
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Fig.  5:  Graph  of  temperature  drop  in  a  thermoelement 
time,  at  various  current  densities  v  (6o=0.6) 


B  W 


Fig.  5:  Continued 
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The  values  of  5^  in  the  case  of  heating  are  lower  than  in  cooling 


at  the  same  values  of  \>,  m  and  Bi*.  For  determination  of  the  value  of 
&l  and  62*  the  data  of  Table  1  can  be  used.  The  difference  Bij-v  should 
be  taken  here  as  the  input  parameter. 


Let  us  consider  further  a  nonequilibrium  thermal  process  in  a 
thermoelement, when  the  temperature  changes  with  time  on  the  hot  side, 
and  account  is  taken  of  heat  transfer  from  the  hot  side  Bi2  and  heat 
capacitance  of  the  hot  junctions  n2-  1°  this  case,  the  boundary  conditions 

for  Eq.  (la)  have  the  form  (1(^=0) : 


(30 

Ol. 


L. -b+*w+ ■*'<»- '■‘JL 

3r  L  “  \*  -  W  -  B,=(5  -  *•>]  l .  • 


(18) 


In  solving  the  problem  by  the  operator  method,  let  us  use  the  Laplace 
transform  for  the  temperature  of  the  cold  junction: 

\rP  | p  (v0„  —  .;2c)  -F  (-.  —  r.i2)  .2i  (ch  Y p  — j)  — 

•  r.?/>)  sh  Y  P _ 

(19) 


0(n)^r  —  ~  £*_ f!l  ]_£  —  gj 0  (•<  —  Bi;  • 

P  '  p*  {  |(v  '•  Blj  c,iP)(v  —  Bi- —  y.?p)  —  p\  sliJ'rp" 


Yp  |Bi,  -r  Bi>  ••  (v  r,j)  p ]  :h  Yp  j 


Expression  (19)  contains  a  pole  at  point  p=0  and  poles  corresponding 
to  solutions  of  equation 


!(-'  +  Bi,  -f  r.i P)(>  -  Bio  -  ■  w)  ~  P\  sh  ~ 

-  V~p  JBi,  H-  Bi,  -j- (r,.  -|-  r^p\  ch  \  p  r-. .  0. 


(20) 
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If  the  thermoelement  parameters  are  such  that  Eq.  (20)  has  only 
purely  imaginary  radicals,  converting  to  the  provisional  domain,  we 
obtain  the  relationship 

A/pn\._  ft  (Bi;  —  y)  (ylip  —  0.5v-)  — 

J (r°)  —  Jo  ( Blj  j.  Bio)  -i-(Bt,  -i-vldj-v)  + 

■  0 V 

1CkV'k  Si,'J5K  ~  K  COr'\  '  C«-.SlS,aJ  S 


Where  CKl  -(v-Bia)v*-^K--w3): 

C,:  --  V-  -f  v0o (V  -  Bio)  -f  o;r  ;/J0; 

CK>  ==  2  -j-  (v  -i-  Bq  -  V},5i)  (2 •/*  +  1 )  -  (v  -  Bio  -j-  v(2o;)  (2rn  -!- 1 ) : 
CKt  ~  Bij  -f-  Bi2  —  (<),  -f-  v;2): 

Cs,  -  o;-  +  4  5'k (V)-.  +  r,2)  -  ±-  (Bi,  -  Bio)  + 

-f  (v  -j-  Bi,  -  vj,0*)  (-1  —  Bi2-f ''JoO*); 

6  is  the  positive  root  of  the  equation 


i{r  i  -  £ _ g|,t.Bla-»(/.1-"ta)  . 

^  (M*5  —  V—  Bi,)  (  BI,  —  r,  O'.-  —  v) 


It  is  interesting  to  compare  the  transition  processes  defined  by 
formulas  (21)  and  (22)  with  curves  which  satisfy  relations  (5)  and  (6) 
obtained  earlier.  It  should  be  noted  that,  as  or  02**®,  the  solution 

of  (21)  and  Eq.  (22)  approach  the  limits  of  (5)  and  (6) . 

The  difference  in  the  steady-state  values  of  the  temperatures  found 
by  (21)  and  (5)  is  equal  to 

?rt  __ _ 0.5v- -j  -<Ou(Bi, -=- •«) _ 

'  ( 1  -I-  v  Bi,)  I  Bi,  Bi;  (v  B:,)(Bu  —  v)) 

and  is  always  positive  when  the  inequality 

v (Bi,  -  Bi,)  —  (Bi,  H-  Bi2  +  Bi,Bi2) <0  (23) 

is  observed. 


It  can  be  shown  that  condition  (23)  is  necessary  for  the  limiting 
solution  of  Eq.  (la)  in  limiting  conditions  (18).  Consequently,  in  accom¬ 
plishing  limitations  on  the  current  feed  (23),  the  finiteness  of  the 
magnitude  of  heat  transfer  leads  to  a  reduction  in  ti.e  effectiveness  of 
thermobatteries  in  the  steady-state  mode  of  operations  in  comparison  with 
the  case  when  812-*®. 


To  evaluate  the  effect  of  values  of  Bi.2  and  02  on  'the  rate  of  the 
transition  processes,  the  values  of  the  roots  of  Eq.  (22)  and  (5)  should 
be  compared.  At  sufficiently  high  values  of  Bi2  and  n2  ft  turns  out  that 
increase  in  value  of  heat  transfer  Bi2  from  the  hot  junctions  leads  to  an 
increase  in  the  radical  6^,  that  is,  a  decrease  in  the  time  for  the  transi¬ 
tion  processes.  An  increase  in  the  magnitude  of  the  thermocapacitance 
n2  effects  the  magnitude  of  6^  similarly  under  the  condition 

y_v  +  Bii 


1/  v  4-  Bii 

»  ’ll 


it 

<T 


and  has  the  opposite  effect  when 


’ll 


it 

>T* 


In  the  general  case,  the  expression  for  a  nonequilibrium  temperature 
trend  with  a  multistage  current  can  be  presented  only  in  the  form  of  a 
series;  this  circumstance  strongly  hampers  its  analysis.  In  particular, 
the  question  of  the  relationship  of  the  extreme  value  of  the  temperature 
drop  in  a  nonequi librium  cooling  process  and  the  current  supply  strength, 
and  the  relationship  between  this  value  and  the  maximum  temperature  drop 
in  a  steady-state  mode  of  operation  remains  unexplained. 

In  order  to  obtain  the  expression  A9  as  a  function  of  time,  it  is 
convenient  to  use  the  model  of  a  thermoelement  with  an  infinite  junction 
length.  Here,  it  should  be  taken  into  account  that,  in  a  sufficiently  short 
duration  of  the  cooling  process  after  switching  on  the  curient,  the  tempera¬ 
ture  field  at  the  hot  junction  still  has  not  begun  to  affect  the  magnitude 
of  the  temperature  of  the  cold  junction.  As  long  as  the  Fourier  number 

clt 

Fo=  is  much  less  than  unity,  the  temperature  trend  on  the  cold  junction 

does  not  depend  on  the  length  of  the  thermoelement  junction.  Therefore, 
it  can  be  considered  infinite.  In  this  case,  the  boundary  conditions  of 
the -heat  conductivity  equation  (l)-(6)  are: 
at  x=0 


OT 


or 


ot 


dT 

tlx 


^0. 


(24) 


(25) 


The  boundary  condition  or.  the  cold  junction  takes  into  account  the 
Peltier  effect,  heat  transfer  to  the  surrounding  medium,  heat  capacity  of 
the  cooling  body,  the  Joule  effect,  separable  at  the  contacts,  and  the 
output  of  the  constant  heat  source. 

First  of  all  let  us  consider  the  case  when  the  cold  junction  is 
isolated  from  the  surrounding  medium  [10];  then  at  x=0 

1  47-  — 
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Let  us  reduce  Eq.  (1)  to  dimensionless  form,  during  which  we  csa  the 
relation  j  __  b  as  the  equivalent  length.  Here, 

<?/ 


d9 


#9  ■  ,  . 

&Ll  '  0(Fo.) 


(26) 

(27) 


<)9 


when  /-  ,>/  tj< 

,  <50  ~ 

when  /— co  -^-—0. 


(28) 


Here 


X= 


jc  _  e/.r  _  of  _  e^at 


d a 


d2 

3 


/  2 


Eq.  (26)  is  easily  solved  by  the  operator  method.  As  a  result  of 
temperature  change  of  the  cold  junction  with  time,  the  following  is 
obtained: 


A9(Fo)=(l  +  0G)(1  —  exp Fo,erfc  V FOj  )  2  ]/-'§*■  . 


(29) 


In  this  manner,  temperature  change  of  the  cold  junction  depends  only 
on  parameter  Fc3,  that  is,  on  the  product  j2t.  Expression  (29)  has  a 
maximum  at  a  specific  value  of  the  argument  Fo3=Fo|  determined  by 
the  equation 


}/ -  Foj  «xp  Foj  erfc  \/  f0;  =  — 


(30) 


Thus,  at  9n=0.6  the  value  of  Foi=0.081  and  AO  =0.0895. 
*  u  3  max 


An  increase 


in  current  j  causes  a  corresponding  decrease  in  the  time  of  arival  of  the 
maximum,  and  vice  versa.  The  value  of  the  maximum  A6  does  not  depend  on 
the  current.  The  greatest  temperature  drop  on  the  element  in  the  steady- 
state  mode  of  operation  is 

“l+0o  --  V 1  +  20o  ,  when  %  --  0,6  A0,1>kC  —  0,1 17. 


The  maximum  temperature  drop  in  the  nonequilibrium  mode  of  operation 
is  less  than  the  maximum  temperature  drop  under  steady-state  conditions. 
Thus,  at  60=0.6,  the  ratio  between  these  values  equals  0.766.  The  data 
shown  are  correct  for  a  thermoelement  with  finite  junction  length  d  if 


the  Fourier  number 


at 

d2 


«  1.  The  maximum  temperature  drop  on  cold 
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junctions  for  a  semi-infinite  element  approaches  at  Fo3=Fo* 
that  is,at-£— (-^)*Fo;. 


e'-par 

7? 


In  this  manner,  condition  « 1  leads  to  the  inequality  v>)'Fo"\ 

As  was  shown  above  [Table  2,  formula  (15)],  the  minimum  current  value  v0 
at  which  an  extreme  appears  amounts  to  0.791,  and  }'Fo*=0,28. 


Beginning  at  the  value  v=v0,  the  values  of  the  temperature  drop  at  the 
extreme  point,  considered  for  models  with  junctions  of  finite  and  semi- 
infinite  lengths,  practically  coincide.  The  trend  of  temperature  change 
A6CTai4  vs*  v  £curve  and  the  values  of  A0  at  the  extreme  point  of  the 

function  A6(Fo)  (curve  2)  for  a  thermoelement  with  a  finite  junction 
length  are  shown  in  Fig.  6.  At  v=vq,  the  maximum  appears  as  t-*»;  with 
increase  in  v,  the  value  of  Fo*  decreases,  and  as  Fo*-*-0. 


Fig.  6:  Temperature  drop  on  a  thermoelememt  with 
thermally  isolated  cold  junction  vs. 
current  density  (6o=0.6) 

Key:  1.  In  the  steady-state  mode  of  operation 

2.  Maximum  in  time 


For  the  general  case,  when  the  limiting  condition  on  the  cold  junction 
is  determined  by  relationship  (24) ,  it  also  is  not  difficult  to  obtain  an 
analytical  expression,  giving  A0  vs.  time.  However,  here,  A0  already  is 
not  only  a  function  of  dimensionless  parameters  0<j  and  Fo,  but  of  the 
dimensionless  complexes 
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The  effect  of  thermoresistances  on  a  junction  decreases  with  increase 
in  current,  and  the  effect  of  heat  capacity  and  electrical  resistance 
increases. 

Fig.  7  presents  A0  (curve  1)  and  A0  at  the  maximum  point  (curve  2) 

CT  014 

vs.  current  v,  for  the  case  Bii=5  (m=0,  £..=0,  K  =0,  6o=0.6).  The  minimum 

current  value  v0  at  which  a  maximum  appears  amounts  to  0.37,  and  here 
vq<v®.  In  the  case  being  considered,  the  temperature  drop  at  the  maximum 

point  is  a  kittle  more  than  at  the  steady-state.  As  the  value  of 

A0  coincides  with  the  temperature  drop  at  Bi=0.  Curve  2  is  plotted  from 
data  calculated  on  a  digital  computer;  the  results  of  the  calculation 
according  t<^  the  formula  for  a  thermoelement  with  a  semi-infinite  junction 
gives  practically  identical  results,  starting  with  v**l.  The  presence  of 
thermocapacitance  and  contact  resistance,  as  has  already  been  demonstrated, 
comes  down  to  the  fact  that  the  value  of  A6  at  the  maximum  point  falls  with 
increase  in  current.  In  the  case  of  the  presence  of  resistance,  the  extreme 
occurs  only  at  sufficiently  small  currents  when  £K<90  or  • pA9q  , 

°r 

when  the  magnitude  of  the  Peltier  heat  exceeds  the  Joule  effect  on  the 
cold  Junction  contacts  at  the  initial  moment  of  time. 


ofl 


Fig.  7:  Temperature  drop  on  a  thermoelement  vs.  current 

density,  taking  heat  transfer  on  the  cold  junction 
into  account  (0q=O.6,  Bii=5) 

Key:  1.  In  the  steady-state  mode  of  operation 

2.  Maximum  in  time 
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